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I. INTRODUCTION 
A. Background and Motivation 

Trapping of neutral antihydrogen was first achieved in 
2010 by the ALPHA collaboration Q] and, by 2011, AL- 
PHA had reported confinement times up to 1000 s [2J. 
Focus is now shifting from proof-of-principlc produc- 
tion and confinement toward precision measurements and 
tests of fundamental physics. 

There are multiple long-term goals motivating antihy- 
drogen research: the first is to search for possible CPT 
violation by examining the spectra of anti-atoms. A first 
step in this direction was taken in 2012, when ALPHA 
measured the frequency of transitions between hyperfine 
levels to a relative precision of 10 -3 [3]. Future work will 
concentrate on high precision measurement of this hyper- 
fine splitting, and of the two-photon IS — > 2S transition. 

A second goal is to search for violations of the weak 
equivalence principle — the equality of the inertial and 
gravitation mass of any object, independent of its com- 
position or structure. Initial experiments with sensitivity 
to gravitational effects of the earth on neutral antimatter 
have been conducted |H [5] , and others are planned [3 U\ ■ 
The ALPHA collaboration has inferred limits on the ratio 
T = — of gravitational mass M to inertial mass m of an- 
tihydrogen by carefully comparing the simulated and ex- 
perimentally determined temporal and spatial distribu- 
tions of antihydrogen annihilations observed during the 
slow turn-off of magnetic trap which confined the anti- 
atoms. Values of T > 110 and F < —65 were rejected [I] 
at a statistical significance level of 95%. In a 2012 publi- 
cation on antihydrogen trapping 5 , the ATRAP collabo- 
ration reported a gravitational measurement that rejects 
T values greater than 200 at a 2 a level. Their methodol- 
ogy was mentioned only briefly, but is apparently based 
on counting annihilation events during radial field shut- 
down in their vertical trap. Understanding the possibili- 
ties and limitations of these techniques on trapped neu- 
tral antimatter provides the motivation for the present 
work. Two other experiments intending to test the weak 
equivalence principle, AEGIS [6] and GBAR [7] also rely 
on the Antiproton Decelerator (AD) at CERN, but will 



use beams of antihydrogen rather than trapped popula- 
tions, and hence their operation is beyond the scope of 
this paper. 

We present here a detailed study of the influence of 
gravity on the nonlinear classical dynamics of trapped 
antihydrogen and, in particular, how features of the non- 
linear dynamics impact gravitational measurement tech- 
niques in vertical traps. Horizontal traps will be dis- 
cussed in more detail elsewhere. Analysis is performed 
in some generality, but specific numerical examples are 
motivated by what we infer are the methodology and, 
roughly, the field geometry used in [5] , as well as the pos- 
sibility of using a vertically-oriented version of the AL- 
PHA trap. Unless we have misunderstood the ATRAP 
methodology, our simulations and analysis show no ef- 
fect of gravity at the levels of sensitivity claimed in their 
measurement, or, indeed, at lower sensitivity. 

B. Dynamical Framework 

ALPHA [8] and ATRAP g] trap neutral antihydro- 
gen in a quasi-static magnetic minimum created by three 
sets of external coils: two mirror coils produce a field B m 
confining the anti-atoms axially, and a multipole coil pro- 
duces a field Bp confining them radially. Both the mir- 
ror and multipole fields will exhibit spatial variation, and 
may also have time dependence. These trapping fields are 
superimposed on a static, uniform background solenoidal 
field B\y = B^z (a legacy of the charged-species trapping 
preceding antihydrogen production), which reduces the 
effective trap depth but also, felicitously, tends to sup- 
press non-adiabatic spin flips in the neutral anti-atoms 
near the field minimum. The total magnetic field at a 
position r and time t is then given by the vector sum 

B{r,t) = B h + B p {r,t)+B m {r,t). (1) 

The orientation of JE?b, which here establishes the lon- 
gitudinal (z) axis, will be assumed to be either verti- 
cal or horizontal, i.e., parallel to or perpendicular to the 
Earth's local gravitational acceleration g. 

The present analysis presupposes a number of other 
simplifying assumptions. A semiclassical approach de- 
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scribes the internal states of the anti-atom quantum me- 
chanically while treating the center-of-mass (COM) de- 
grees of freedom classically. While we allow for the pos- 
sibility of matter/antimatter asymmetry in gravitational 
interactions, such that the anti-atom may have an ef- 
fective gravitational mass M different from its inertial 
mass m, we shall here presume, consistent with CPT 
invariance, that each antihydrogen is precisely electri- 
cally neutral, so experiences no net external Coulomb 
nor COM Lorentz forces, but can experience a Zeeman 
force due to a non-zero expectation value for its mag- 
netic moment, whose magnitude is identical to that of 
an ordinary hydrogen atom in the corresponding inter- 
nal quantum state. 

Trapping times are sufficiently long so that we may 
confine attention to anti-atoms which have relaxed to the 
orbital ground state |10j . in which the anti-atom's mag- 
netic moment is dominated by the positron spin, and 
to low-field-seeking spin states, such that the anti-atoms 
can be trapped near a magnetic minimum. (Those anti- 
atoms in a high-field seeking state quickly hit the trap 
wall and are thus not considered here). The effects of 
the magnetic field on the internal state are small since 

H B B = \hVL e < a 2 mc 2 , (2) 

where a = e 2 /hc (in Gaussian units) is the fine structure 
constant, h = h/2n is the reduced Planck's constant, c 
the speed of light in vacuo, m is the rest mass of the 
positron and e is the magnitude of its electric charge, 
while /ib = eh/2mc is the Bohr magneton, and Q e — 
Q e (r,t) — eB/mc is the local positron gryofrequency, 
where B = B(r,t) = \B(r,t)\ is the magnitude of the 
total magnetic field at position r and time t. 

Characteristic antihydrogen translational tempera- 
tures T are such that translational motion remains en- 
tirely non-relativistic: 

fceT <C tyic , 

where here fee is Boltzmann's constant. The tempera- 
ture should also be sufficiently low so that as the anti- 
atom translates, the changes experienced in local mag- 
netic field strength remain adiabatic with respect to spin 
dynamics: 

V / ¥^ M «^- (3) 

for spatial positions r accessible to trapped anti-atoms. 
A related asumption is that the characteristic radial 
bounce frequency uj r and longitudinal bounce frequency 
to z are both small compared to J7 e . Under these assump- 
tions, the magnetic moment adiabatically tracks the di- 
rection of the field, and the classical COM dynamics is 
governed by the Hamiltonian 

H = H(r,p,t) = ^p 2 + [i B B(r,t) + Mg-r, (4) 

where p = \p\, p = mr is the momentum of the anti- 
hydrogen. The + sign in front of /x^ is appropriate for 



anti-atoms that are in a low-field seeking hyperfine state 
— those that can be stably trapped in a local minimum 
of B(r,t). 

The study of antimatter gravitational forces requires 
examining the antihydrogen dynamics for various as- 
sumed values of the ratio T = M/m and comparing these 
results with experimental observation. The gravitational 
force modifies the dynamics in ways that depend on trap 
orientation, on field geometry, on initial conditions, and 
on the time profile of the trap turn-off. In ALPHA, an 
octupole field provides transverse confinement, and the 
trap axis is horizontal, perpendicular to g [S]. ATRAP 
instead employs a quadrupole for transverse confinement, 
while the trap axis is vertical, parallel to g [9]. 

In cylindrical coordinates (r, <f>, z) oriented such that 
the solcnoidal field points along z, the squared- 
magnitude of the total field ([I]) is 

B 2 = B 2 (r, <j>, z, t) = B 2 r (r, cp, z, t) + B 2 (r, 0, z, t) 

+ [B m Z (r,z,t) + B pz (r, 0, z, t) + B h } 2 , (5) 

where B r and B^ are the radial and azimuthal field com- 
ponents, which arise from both the multipole and the 
mirror coils. Clearly, the general form of the trapping 
potential \ibB will depend non-trivially on (r, </>, z) and 
possibly on t if field turn-off is modeled. However, since 
the decay of the magnetic fields is very slow, unless oth- 
erwise noted, our dynamical analysis will be performed 
presuming a frozen value of the Hamiltonian, and the ex- 
plicit t dependence in B = B{r) or therefore H = H(r,p) 
will generally be suppressed in the notation. 



C. Basic Dynamic Considerations: Regular versus 
Stochastic Trajectories 

We can gain some basic understanding of the dy- 
namics if we temporarily ignore the radial component 
of the mirror field and any end effects from the multi- 
pole. Under these simplifications, an ovder-i multipole 
field yields B 2 + B 2 oc r 2e ~~ 2 , with no ^-dependence 
in field magnitude, i.e., B = B(r,z). For a vertical 
trap, the total (magnetic and gravitational) potential 
U(r) = hbB + Mg ■ r is then axially symmetric and 
separable, i.e., U(r) — U r (r) + U z (z), and the longitu- 
dinal antihydrogen motion along z is uncoupled from the 
transverse motion in the (r, (f>) plane. The trajectories 
are regular and fully determined by integrating two one- 
degree-of- freedom Hamiltonian systems, namely 

H\\(z, Pz ) = ^p 2 z + U z (z), (6a) 
H ± (r,p r ; P4> ) = ^pl + j^jpl + U r (r), (6b) 

where p r and p z are, respectively, the radial and longitu- 
dinal components of the momentum, and the azimuthal 
component p^ represents the angular momentum along 
z. 
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In this situation, there are three dynamical invariants: 
the perpendicular energy H± = E±, the parallel energy 
H\\ = , and the angular momentum — mrv^, where 
v 4> = iv x ~ iy)/ r is the azimuthal velocity. Any one 
trajectory with total energy E = Ei\ + E± will not be er- 
godic and will not explore the entire energy hypersurface 
H(r,p) — E. Consider the consequences of very slowly 
lowering the transverse confining field, i.e., U r — > as 
t — >• oo. With uncoupled motion there will be no correla- 
tion between the axial dynamics, where gravity acts, and 
the transverse dynamics; the motions are uncoupled and 
non-ergodic. Consider an anti-atom that has a kinetic en- 
ergy below the level of transverse potential barrier, but 
above the axial potential barrier. Such an anti-atom is 
confined transversely, but may escape the trap axially. 
Nonetheless, it will remain confined if enough of its en- 
ergy is tied up in transverse motion; because there is no 
coupling, it would never come to possess sufficient energy 
to overcome the axial barrier. 



II. ANALYTICAL DESCRIPTION OF SINGLE 
ANTI-ATOM MOTION 

Detailed analysis of single anti-atom dynamics in a 
magnetostatic trap is crucial for understanding antihy- 
drogen losses, laser cooling of trapped antihydrogen, and 
limitations of different approaches to measuring the grav- 
itational mass of antihydrogen. While there is no general 
solution for the full three-dimensional anti-atom trajec- 
tory in arbitrary fields, the analysis can be considerably 
simplified and some insight provided by the case of a trap 
with nearly-separable confining potentials. 

To begin, we apply canonical Hamiltonian perturba- 
tion theory |13j to analyze single anti-atom motion and 
then use obtained results in Sec. |III C] to compare numer- 
ical simulations with analytical predictions. 

A. Perturbation theory for weakly coupled motion 



In more realistic geometries, some amount of coupling 
will be caused by the radial components of the mirror 
field, by the octupole end effects, and by higher-order 
multipole contributions. If the coupling were sufficient to 
make the motion fully crgodic, then an antihydrogen with 
total energy E exceeding the lowest of the axial trapping 
potentials U m - m would eventually escape. Knowing U m - ln 
would then allow one to obtain a bound on M by vary- 
ing only the radially confining potential. Such ergodicity 
may have been implicitly assumed in the gravitational 
discussion in [5]. 



Since the dynamics in a realistic magnetic field formed 
by mirror and multipole coils are not fully integrable, nor 
expected to be fully ergodic, numerical simulation may be 
required. For small coupling, large regions of phase space 
should remain integrable. The KAM theorem [TT] sug- 
gests that the majority of resonant tori will survive suf- 
ficiently small perturbations and the corresponding tra- 
jectories remain quasiperiodic. In this case, many anti- 
atoms would remain trapped even when, ostensibly, they 
appear to have sufficient energy to escape axially. 



The remainder of this paper is organized as follows. In 
Sec. |TlJ a perturbation theory is used to study the influ- 
ence of coupling on the anti-atom dynamics. A discussion 
of numerical issues and detailed simulation results for a 
vertical trap with a field profile similar to that of the 
ATRAP experiment are presented in Sec. Ill These re- 



sults indicate weak coupling between the transverse and 
longitudinal dynamics, which may be typical for other ex- 
isting atom traps [H] as well. An alternative approach to 
measuring antihydrogen gravitational mass in a vertical 
trap, which involves turning off the mirror fields (which 
does not require anti-atom ergodicity), is also discussed. 
Our conclusions are given in Sec. |IV| 



Consider antihydrogen motion in an almost separable 
potential U(r,(f>,z), i.e., U(r,(j),z) = U r (r) + U z (z) + 
8U(r, (j>, z), with SU much smaller in magnitude than 
U r (r) + U z (z\. After rewriting the original (frozen) 
Hamiltonian El as 



ff(r,0,z,p r ,p0,p z ) 



Pr 

2m 



2mr 2 

+ U z (z) + U r (r)+6U(r, tf>,z) 



2m 

= H + SU, (7) 



the term SU can be considered as a perturbation to the 
integrable system with integrable Hamiltonian H . For 



trap designs like those to be discussed in Sec. IIIB the 
magnitude of the perturbation SU may be comparable in 
relative magnitude to the trapping well depth (reaching 
0(1/3) for the vertical quadrupole trap). However, at 
least in our examples, this maximal value for SU is ac- 
cessible to only those anti-atoms that are weakly trapped 
in both radial and axial directions, so the perturbation 
theory should provide some insight into more typical tra- 
jectories. 

The first step in applying the perturbation theory to 
Eq. ([7 ) is to find the action-angle variables of the unper- 
turbed Hamiltonian Hq. The axial motion is uncoupled 
from the transverse oscillations |13j : 



i 

2^ 



p z dz, 



(8) 



where I z is the axial action, and the integration is per- 
formed over a closed trajectory solving pl/2m + U z (z) = 
H\\ . The frequency of the axial oscillations is lj z = ip z = 
dHu/dI z , where ip z is the angle variable canonically con- 
jugate to I z . Assuming that the magnetic field profile 
U z is almost quadratic in some vicinity of z = 0, u> z {I z ) 
is nearly constant for small I z (when the anti-atom os- 
cillates near the center of the trap). However, when I z 
becomes so large that the anti-atom trajectory passes 
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FIG. 1. Phase space portrait of dynamical system governed 
by Ha — p 2 z /2m + U z (z) showing two separatrices passing 
through the lower Si (bold solid) and the higher S2 (dot- 
ted) axial potential barriers (top); and a corresponding ax- 
ial potential profile U z (z) (bottom). The axial coordinate 
z is normalized to the device half-length L, p z is normal- 
ized to {2mma.x[U z (z) — C/ z (0)]} 1//2 , and U z is normalized to 

wax[U g (z) - U x (0)]. 



near one of the mirror coils, the frequency Li z (l z ) de- 
creases, vanishing eventually at I z = I*, where the anti- 
atom turning point reaches a local maximum of U z (z). 
This trajectory corresponds to a separatrix in the (p z , z) 
phase space (Fig.[l]). 

The radial action I r can be found similarly. First 
consider a canonical transformation of H± = p 2 ./2m + 



U r (r) effected by the generating function: 



#(r, </>, I r , P4) = / Pr(r; I r , P<f,) dr + Prf, (9) 



where p r (r; I r , Pq) solves equations of motion for p^ = P^ 
and Hx = H±(I r , P$), with Ir{Hj_,p^) given by 



Ir(H A 



2n 



p r (r;iJj_,p0)dr. (10) 



After this transformation, H± becomes a function of the 
new actions I r and and is independent of the new 
27r-periodic angles ip r and ip^, = <p + A, where 



(11) 



The canonical angles are generally defined up to an over- 
all constant. In the following, we choose ip r = when the 



anti-atom is closest to the device axis and ip z = when 
the z coordinate reaches its maximum, and A(r = 0) = 0. 

It is generally difficult to obtain analytical expressions 
for the frequencies ui r = tp r = dH±/dI r and = tp^ — 
dH±/dP c f > . Their values are related at P^ — 0, when 
the anti-atom velocity has a vanishing azimuthal compo- 
nent. Introducing the full period of transverse oscillation 
T±(I r ), one can see that <p has a period T_l, while ip r has a 
period T±/2. Then, recalling that <p = ip^,—A(ip r ; I r , P4,), 
one obtains ip^{t + T±/2) - ip^t) = <j>(t + T ± /2) - <p{t) 
and, therefore, lo^T^/2 = ir, or lu^ = cu r /2. 



B. Axisymmetric perturbation 

First, consider the case of a purely axisymmetric per- 
turbation 8U{r,z). The Hamiltonian written in action- 
angle variables is 



H(I,iP) = H (I r ,I z ) 

OO OO 

+ E E SU ktl (I r ,I z )e ik ^ 



-ilip z 



(12) 



where 8Uk,i are the radial/azimuthal Fourier components 
of 5U{xp\I). The only resonances are, therefore, of the 
form kuj r + Iuj z = 0. This can be rewritten as: 



dz 



y/H^-U„(z) 



dr 



(13) 



E-B«-U r (r)- 



2mr 2 



where z\, Z2, ri, and ri are the longitudinal and radial 
turning points, and Hu is a function of I z . Note that for 
the pure quadrupole field with U r ~ r 2 , the right-hand 
side of Eq. (13) is independent of E, I z , and P^, while 
the left hand side is a function of I z only. 

Consider a long trap with the radius R w much smaller 
than the longitudinal half-length Z. The ratio U) r /oj z 
scales as Z / R w and the axial scale of the perturbation 
SU will be on the order of the characteristic coil ra- 
dius R c . Since the perturbation maximum is reached at 
R « R w near the mirrors, the resonance harmonics SUk,i 
and the corresponding resonance widths SI ~ y/\SUk,i\ 
grow with increasing anti-atom energy. But even for 
the highest-energy anti-atoms, 511^ 1 is roughly propor- 
tional to (R c /L) max[5£7] exp(— k\R c /R w \) and is small. 
In this case, the radial anti-atom oscillations are adi- 
abatic [HI [15]. Assuming that most neighboring res- 
onances do not overlap, the system dynamics within 
resonance islands is expected to be regular, becoming 
stochastic in small vicinities of the island separatrices 
only. However, since oj z {I z ) vanishes at the critical point 
I* (see Sec. II A ), there will be an area in the phase space 
where resonances accumulate and overlap [TB] , thus form- 
ing a stochastic layer in a vicinity of the separatrix at 

iz = r z mm- 

The systems with axisymmetric perturbations 8U[r, z) 
possess another non-generic property which does not 
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survive once azimuthal angle dependence is introduced. 
Namely, for small SU, all antihyrogen trajectories, even 
stochastic ones, are localized and cannot change their 
energy by more than a certain finite amount. With any 
amount of angular dependence, this is no longer the case 
and, in fact, some trajectories in certain time-dependent 
two-dimensional dynamical systems are known |17j to 
"diffuse" indefinitely albeit slowly, reaching any cho- 
sen value of action at some sufficiently late moment of 
time. Known as Arnold Diffusion, this phenomenon will 
be at least partially responsible for slow anti-atom loss 
from the stationary quadrupole and octupole traps (see 
Sec. IIIIC lb. 



C. Non-axisymmetric perturbation 

1. Resonances 

Consider next an angle-dependent perturbation of the 
form SU(r, <fr, z) — 5V{r, z) cosncj) for some fixed integer 
n. The Hamiltonian is now: 



H(I,i/>) = H Q (I r ,I s ) 

k,l 



-incp\ ikil> r -\-ilip z 



(14) 



where 8Vk,i is the angular Fourier component of 
dV(if>; I). After substituting (j> = i/j^ — A(ip r ; I r , P<f,), 



H = H {I r ,I z ) 

+ Y.{ 5w ^ I ^ I *)^ r+l ^ z+ 



k.i 



+ W M ,_ n (J r , I z ) e^+^'-^j , (15) 

where SWk,i, n is calculated given all the SVk,i as well as 
A(?/v). Since the angle-dependent harmonic is fixed, the 
resonance condition becomes: 



Qk,l,±n(I) = kuj r + lu z ± = 0. 



(16) 



Since all frequencies in Eq. (16) are functions of I r , I. 



and Pa,, one can find the resonances in action space. Fix- 
ing the total energy H {I r , I z , P<p) — E, the resonance 
curves can, for example, be plotted in the (I z ,P<p) coor- 
dinates. Let 1Z„ be the set of such curves in the (I z , P^) 
space corresponding to Qk,i,± n — for some k, I G Z. 
Such a plot is shown in Fig. [2] for the quadrupole trap 
design discussed in more detail in Sec. Ill C Notice that 
the resonances ku> r + lu z = in Fig!~[2[ are character- 
ized by nearly constant values of I z , due to the fact that 
both sides of Eq. ( 13 ) are independent of P^ for the pure 
quadrupole field. 
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FIG. 2. (Color online) Resonances kui r + Ilj z + suj^ — 
in the {I z ,P<t>) space plotted for E — 390 mK in a vertical 
quadrupole trap: jsj = (solid red), |s| = 1 (dash-dotted 
blue), |s| = 2 (dashed magenta), |s| = 4 (dotted green). Only 
resonances with \k\ < 10, \l\ < 12 and sj < 12 are shown. 
The action I z is normalized to the action of the separatrix 
trajectory corresponding to I z = 1, and = P r /,/(mR w v), 
where v — \/2E /m and R w is the wall radius. Note that for 
a perturbation with n — 4, all sj — 1, sj — 2 and jsj = 4 
resonances affect anti-atom dynamics. 



tant resonances. For sufficiently small 6V, the charac- 
teristic width of the resonance Qk,i,s(Ir, Izi P<j>) = is 
defined by the amplitude of resonant oscillations AI r w 
|/s|AJ, AI Z w |Z|A7 and AP w |s|AJ, where [18] 



AI 



\6Wj 



k.i.. 



\%H \ 



(17) 



and d t = k d/dl r + ld/dl z + s d/ dP^. 

Although, for a wide class of smooth functions, the 
widths AI are expected to decrease exponentially with 
\k\, \l\, and |s|, the calculation of the exact value of 
SWk^i^s is generally quite complex. However, it can be 
simplified for the Qk,i.o = and Qh,o.-2h = resonances. 
Indeed, 8Wk i a f° r n = is given by: 



SWk 



L0 



\5Vu 



k.i- 



(18) 



On the other hand, recalling that Qh.o,-2h vanishes at 
Ps = 0, one obtains for n = 2h: 



2. Resonance widths 

Characterizing anti-atom dynamics in phase space re- 
quires a knowledge of locations and widths of all impor- 



2tt 2tt 



k=-c 



SV k ,o cos(2^) 



^ikifj r — ihip r -\-2ihtp ( i ) 



dtprdip^. (19) 
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After substituting cf> = tp^ — A, this becomes: 



1 



^ fSV k , 



SWh,0,-2h = 7T / i — r~ 
Z7T » — ' J i 

k— — oo q 

• ^ihA+ik^-ih^r d ^ r = sVa fi /2, (20) 

where we have used the fact that 2A = ijj r (mod 2ir) for 
— 0. In the following section, we use Eqs. (161, (181, 
and (20) to find resonances, estimate their widths, and 



reach some qualitative conclusions about antihydrogen 
dynamics in the trap. 



III. NUMERICAL SIMULATIONS 

In this section, numerical simulations of single anti- 
atom motion, aimed at assessing the ergodicity of anti- 
atom trajectories and studying the feasibility of gravita- 
tional measurement techniques, are discussed. 



A. Computational framework 

Antihydrogen dynamics were simulated using both 
standard Runge-Kutta and fourth-order symplectic 
schemes to integrate the COM equations of motion r = 
+dH/dp and p — —dH/dr, where the frozen Hamilto- 
nian H is given by Eq. (El. The magnetic field profile 
B(r) was calculated from the presumed configurations of 
magnetic coils at some fixed reference time. Since deter- 
mining the magnetic field using the Biot-Savart law or 
series expansions for each anti-atom at each moment of 
time would be quite computationally expensive, we pre- 
calculated Bona fixed lattice and then interpolated B 
at instantaneous anti-atom positions. Given the repre- 
sentation B(r) = B2n(r, z) cos(2n0 + 2n ), in the 
configurations of interest, harmonics B n with n > 2 for 
quadrupole traps and n > 4 for octupole traps can be ne- 
glected. Therefore, instead of storing a three-dimensional 
array of B values, we calculated Bo(r,z) and B2\i(r,z) 
on a two-dimensional lattice. The angular harmonics B 
and £?2\4 were calculated using a fast Fourier transform of 
B(r) on a ring (r, <pi, z), where 4>i = 2ni/N with N = 64. 

Using only bilinear interpolation to find Bq and B 2 \4 at 
some intermediate point would be undesirable, since the 
force acting on each anti-atom is proportional to VB, 
which would then be a discontinuous function causing 
noise and large numerical errors in antihydrogen trajecto- 
ries. Instead, we used a bicubic interpolation 19 , which 
produced a C 1 -smooth approximation of B(r). 



B. Trap geometries 

Three device designs were considered: (a) a vertically- 
oriented quadrupole trap with parameters similar to 




FIG. 3. Set of coils used in our numerical simulations: 
quadrupole coils (solid) with current directions shown with 
arrows and mirror coils (dashed). The solenoid creating a 
constant background field is not shown. 



those of the ATRAP experiment, in which ergodicity 
and feasibility of gravitational measurements via lower- 
ing of the radial confining potential were studied; (b) a 
vertically-oriented octupole trap otherwise similar to the 
ALPHA apparatus, which we used to analyze alterna- 
tive approaches to antihydrogen gravitational mass mea- 
surements; and (c) a horizontally-oriented octupole trap 
similar to the actual ALPHA apparatus (to be discussed 
elsewhere) . 

In all device designs, the background magnetic field 
Bb w as directed along z, with magnitude equal to IT. 
In a device design motivated by ATRAP [5j [9] (but not 
an exact model), two mirror coils of radius R p — 10.4cm 
were placed at \z\ = Z = 10 cm (Fig. [3]). The total cur- 
rent of approximately 265 kA flowing through each mir- 
ror coil increased the magnetic field at the trap center 
to 2.2 T, while creating a 375 mK axial well depth. The 
quadrupole coil was modeled as a combination of 4 rect- 
angular loops with longer sides of length 2Z = 20 cm 
directed along z and shorter sides of length Ri sa 6 cm 
directed along either x or y (Fig. [3]) . Each loop coil lo- 
cated at |x| = Ri or \y\ = Ri carried the total current 
of approximately 360 kA. As a result, a 375 mK radial 
well was also created. The trap walls, on which the an- 
tihydrogen are assumed to immediately annihilate, were 
chosen to be at \z\ = 1.2 Z and at r = R w = 1.8 cm. In 
a realistic trap, there are no actual walls at \z\ = 1.2 Z, 
but all anti-atoms reaching this location will never return 
to the trapping volume and will annihilate shortly there- 
after. Performing the angular Fourier decomposition of 
B inside this volume, one obtains \B±(r, z)\ < 60 G, while 
B a is between 2 T and 3 T, and \B 2 \ < 0.25 T. Neglecting 
octupole and higher-order angular harmonics is therefore 
justified for this trap. 

In a trap design based on that of ALPHA, the mirror 
coils located at \z\ = Z = 13.7cm created a 670 mK 
axial well depth for antihydrogen. The octupole coil 
was modeled as a combination of 8 rectangular loops 
with longer sides of length 2Z, located at a distance 
Ri = 2.3 cm from the device axis and connected by 
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shorter sides of length 2Ri tan(7r/8). The magnetic field 
created by the octupole reached 1.5 T on the trap wall 
(R w = 2.22 cm). Simulated anti-atoms were assumed 
to annihilate upon encountering this wall or else when 
reaching \z\ = L = 15 cm. 



C. Vertical trap simulation 

1. Ergodicity of anti-atom trajectories 

As discussed in Sec.|TJ a suggested method of measuring 
the gravitational mass of trapped antihydrogen by low- 
ering the radial well depth of a vertically-oriented trap 
and observing annihilations of escaping anti-atoms [5] 
implicitly assumes the ergodicity of anti-atom trajecto- 
ries. The assumption of ergodic trajectories with energy 
higher than the lowest of the axial trapping potentials, 
U m m ~ 350 mK, was verified numerically by simulat- 
ing anti-atom esc ape fr om the vertical quadrupole trap 
described in Sec. IIIB In our simulations, 6 • 10 3 



In our simulations, 6 • 1CP anti- 
atoms were initialized in the trap center homogeneously 
within a cylinder of radius 1 cm and length 2 cm. The 
ratio T = M/m of the anti-atom gravitational mass M 
to the inertial mass m was chosen to be 200, in accor- 
dance with the limit asserted in [5]. Initial anti-atom 
velocities were distributed isotropically, and their en- 
ergies were chosen randomly and homogeneously from 
a range J7 m ; n < E < 550 mK. The anti-atom phase 
space positions were then numerically evolved in time 
in static fields for 10 3 seconds, during which a typical 
anti-atom performed about 2 • 10 5 axial and more than 
7 • 10 5 transverse oscillations. Every anti-atom encoun- 
tering the device wall was assumed to annihilate imme- 
diately, causing the total number of trapped anti-atoms 
to drop over time. Figure [4] plots the simulated frac- 
tion of anti-atoms remaining in the trap as a function of 
time, f{t) = n(*)/n(0), from t = s to t = 1000 s. Dif- 
ferent numerical integration schemes showed good agree- 
ment, and indicated that after 1000 s more than 25% of 
all anti-atoms remained trapped in the device with only 
2.5% of anti-atoms escaping in the last 999 seconds. That 
is to say, most escaping antihydrogens escape very early 
— about 90% anti-atoms which do escape leave the trap 
within the first 10 ms, which is comparable to a single ax- 
ial bounce time. The fact that in our simulations there 
exist anti-atoms trapped in the system for 1000 s is not 
consistent with the assumption of ergodicity; instead, it 
indicates the existence of bounded regular trajectories. 

It is instructive to see the effect of the angle-dependent 
harmonic in B(r) on the anti-atom escape rate. In 
one of our simulations, we considered an axisymmetric 



potential U(r,z) = [iB q0 ( 



r,z) + 



Mgz identical to the 



quadrupole potential except for an artificially suppressed 
^B q 2{r, z) cos 2<j) term (Fig. El. This potential has the 
same angle-averaged profile as in the quadrupole field, 
but it cannot be physically realized. In this case, similar 
to the previous quadrupole simulation, more than 50% of 




1000 



FIG. 4. Fraction f(t) — n(£)/n(0) of anti-atoms remaining 
trapped in a device as a function of time. Simulations were 
performed for 6000 anti-atoms with energies within a range 
350 mK < E < 550 mK for a trap with a quadrupole coil sim- 
ulated by B(r, </>, z) = B q o(r, z) + B q 2(r, z) cos 20 (blue solid) 
and a trap with an octupole coil simulated by B(r,(f>,z) = 
B o (r, z) + B 4 (r, z) cos 4<j) (red solid) . Antihydrogens escape 
from axisymmetric potentials which, although not realizable 
in multipole traps, share the same angle-averaged profiles with 
quadrupole and octupole traps were also simulated: a trap 
with U (r, z) — fj,B q o + Mgz (blue dashed) and a trap with 
U(r, z) = /jBoo + Mgz (red dashed). The gravitational to in- 
ertial mass ratio T — M /m was equal to 200 in all simulations 
except for one, where anti-atom escape from a quadrupole 
trap with B(r,<f>,z) — B q o{r,z) + B q 2 (r, z) cos 2<f> assuming 
T = 1 was analyzed (blue dotted). 



all anti-atoms escaped within the first 10 ms. However, 
f(t) at later times was much flatter in the axisymmetric 
system, suggesting that the angular resonances may be 
responsible for a slow anti-atom transport in phase space. 
Indeed, such resonances may lead to Arnold diffusion, 
which slowly empties resonance layers, driving resonant 
anti-atoms to the walls. 

The resonance effect of the angular harmonics is even 
more strongly pronounced in fields with higher multi- 
pole perturbations. To observe such effects, we simulated 
anti-atom dynamics in the octupole field by changing the 
total number of loop currents in the vertical quadrupole 
trap described in Sec. |IIIB] from 4 to 8 and also increas- 
ing the current I q by approximately 4 times to create 
a similar radial potential barrier, while also reducing the 
transverse coil size from R to 0.6 R. The survival fraction 
f(t) obtained for the octupole field B = B q q + B Q 4 cos 40 
and the axisymmetric potential U(r,z) = fiB oQ + Mgz 
are shown in Fig. [4] Although the axisymmetric potential 
U(r,z) — /iBoQ + Mgz cannot actually be realized in a 
multipole magnetic trap, simulating anti-atom dynamics 
in it helps to highlight the role of angular perturbations 
in long-time anti-atom dynamics. 
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FIG. 5. Frequency map analysis (FMA) diagram |20| showing D(P^,,I Z ) — In 1 1 — Q,b/£Ia\ for P$(t) in a vertical quadrupole 
trap. Primary frequencies J1a and Qb are calculated numerically over two successive 4 second time intervals A and B. Blue 
and green colors correspond to regular quasiperiodic trajectories. Values of P<$, and I z together with ip r = ip z — ip$ — define 
the initial conditions for the system trajectory with: (a) E = 390 mK and (b) E — 475 mK. 




0.05 0.1 0.15 0.2 0.25 



FIG. 6. A FMA map similar to that in Fig. [5^,, ob- 
tained for E = 390 mK in a system with a trapping poten- 
tial U = Uo + ((72/4) cos 20, where Uq and U2 are calculated 
for the vertical quadrupole trap (Sec. Ill B I . Color palette is 
altered to highlight the parabolic curve (white arrows) con- 
taining strongly nonlinear trajectories. It corresponds to the 
boundary of the resonance island over the stable stationary 
point (Fig. [7k. The straight line at P^ = corresponds to 
trajectories near the saddle point. 



2. Comparison with analytical predictions and frequency 
map analysis 

The majority of anti-atoms trapped for more than 
10 ms in the quadrupole trap configuration considered 
above would appear to exhibit regular trajectories. This 
can be explained qualitatively using the formalism out- 
lined in Secjn] After calculating I r and I z numerically, 



using Eqs. (I 
are obtainec 



and (10), the frequencies ui r , uj^, and w 2 



yy differentiating the unperturbed Hamilto- 



nian expressed as a function of the corresponding actions 
J. Knowing these canonical frequencies, we identify all 
resonances for anti-atoms with fixed energy E and plot 
them in (I z ,P,p) space. Figure [2] shows s uch a plot for a 
vertically-oriented quadrupole trap (Sec. IIIB) and E = 
390 mK, with resonances Qk,i. s = koj r + lu) z + sou = 
and \k\ < 10, |i| < 12, and \s\ < 12. (Without limitation 
on k, I, and s, the entire plot is covered by a dense set 
of curves at this resolution) . The triangular shape of the 
plot is due to the fact that H^_(I ri P^,) = E — Hn(I z ) 
decreases as I z increases. 

Not all resonances shown in Fig. [2] influence the dy- 
namics significantly. For example, consider a trap with 
a perturbation possessing only one angular harmonic, 
i.e., SU(r,(f),z) — V(r, z) cos ncj). For n = 0, correspond- 
ing to an axisymmetric perturbation, all resonances have 
the form ku r + Ilj z = 0. For n = 1, the only reso- 
nances affecting anti-atom motion are those TZf solving 
Qk,i,±i = (shown with dash-dotted lines in Fig. |2b. 
The number of resonances increase with a quadrupole 
field (n = 2). Indeed, every resonance Qk,i,±i = is 
also a resonance for n = 2 since Q2k,2i,±2 = 2Qk,i,±i> 
i.e., IZf C TZ-2 ■ Other resonances Qk,i,±2 = for which 
either fc or / is an odd number, including the resonance 
iu r = 2oj0 at i-0 = (see Sec. II A), are shown in Fig. [2] 
with yellow dashed lines. For the octupole perturbation 
with n = 4, the number of resonances increases even fur- 
ther since, again, TZf C IZf C TZf, i.e., the set of all res- 
onances Qk.i,±A — also includes resonances Qk,i.±i = 
and Qk,i,±2 = 0. This effect may be partially responsible 
for the presence of a larger fraction of anti-atoms with 
stochastic trajectories in the octupole traps (Fig. [4|. 

The fraction of anti-atoms affected by a specific reso- 
nance depends on its width. Using Eqs. (18) and (20), 
the widths of resonances Qk,i.o — and Q/i, 0.-2/1 = can 
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FIG. 7. Dependence of | {P$(t)) \ (T z , Pp) calculated for E — 390 mK in a system with a trapping potential U = C/o+(^W4) cos2<^, 
where Uq and U2 are obtained for the vertical quadrupole trap (Sec. |IIIB" 1: (a) all anti-atoms are initialized with ip r = ip^ = 
ip z — 0, (b) anti-atoms are initialized with ip$ = tt and ip r — ipz — 0. The boundary (dotted line) of the resonance uj t = 2uj$ at 
— calculated using Eq. ( 20 1 contains two branches, for which the corresponding one-dimensional phase portraits are shifted 
by 7r with respect to each other. For one of these branches [bottom curve on figure (a)] the initial condition with ip r = ipj, = 
initializes the system state over the saddle point, hence (P^) is non-zero. For the other branch, this initial condition places the 
system state over the stable stationary point, which makes (Pj,) = within the maximum width of the resonance island. 



be calculated numerically for a vertically-oriented trap 
with parameters similar to those of ATRAP (Fig. [2]). The 
resonance <3i,o,-2 = is then shown to affect a large frac- 
tion of trapped anti-atoms, while Qk,i,a — resonances 
have much smaller widths. Therefore, since there is no 
resonance overlap over a large phase space volume, most 
anti-atom trajectories are expected to be regular. 

The predicted locations of resonances (along with their 
widths) and the associated stochastic layers can be veri- 
fied numerically using a variation of the Frequency Map 
Analysis (FMA) method [5D]. The core idea behind this 
technique is to test system coordinates like Ike l ^ k for 
quasipcriodicity. Treating such a variable as a function 
of time, one can approximate it as a sum of harmonics 
J2k=i Ake lQkt with < \Ak\ and then compare the 

values of Ak and f2fc on different non-intersecting time 
intervals. If the frequencies and amplitudes change con- 
siderably along a single anti-atom trajectory, it can be 
regarded as a sign of stochasticity. Figure [5] shows the 
FMA maps obtained by analyzing P^,(t) for initial an- 
gles ip<p = ip r = ip z = and energies E — 390 mK and 
E = 475 mK. The initial conditions corresponding to 
chaotic anti-atom motion are shown with yellow and red 
in Fig. [5] This figure and other numerical results ob- 
tained for different initial angles suggest that nearly all 
anti-atoms with energy E = 390 mK are characterized 
by quasiperiodic trajectories (shown with blue) rather 
than chaotic motion. Higher energy anti-atoms with 
E = 475 mK, however, are more likely to exhibit stochas- 
tic dynamics. Emergence of stochastic orbits for higher 
energies can be attributed to the fact that many anti- 
atoms can now reach regions near the wall (r = R w ) at 
z « Z, where the angle-dependent perturbation of the 
trapping potential becomes particularly strong. Note, 
however, that emergence of stochasticity does not neces- 



sarily imply rapid anti-atom loss. In fact, the majority 
of anti-atoms with initial states within the bright area in 
Fig. [HJd were shown to stay in the system for at least 100 
seconds. 

Some of the main features of the FMA maps shown 
in Fig. [5] can be related to our analytical predictions. 
Since agreement is better observed for weaker pertur- 
bations 5U (characteristic of traps with smaller radii), 
consider an artificial system with a trapping potential 
U(r,(f>,z) = Uq + (J72/4) cos 2(f), where Uq and U% are 
calculated for a vertical quadrupole trap discussed in 
Sec. |III B| Calculating locations and widths of resonance 
islands, all resonances except for uj r = can be shown 
to affect only a small region of the system phase space. 
On the other hand, the width of the resonance cu r = 2uu 



calculated using Eq. ( 20 1 is sufficiently large to affect al- 
most half of all 390 mK anti-atoms. Interestingly, the 
perturbation harmonic £^1,0,2 corresponding to this res- 
onance and considered as a function of I z for = 
passes through zero at some I z = J*. This means that 
the phase portrait of the resonance island shifts in phase 
by tt after I z goes through 7*. As a result, both when 
I 2 < i*, = or when I z > /*, ips = 7T, the anti-atom 
is initialized over a saddle point and, thus, all such orbits 
are not trapped within the resonance, but lie outside of 
the resonance island. On the other hand, for I z > /*, 
xptj, = or I z < i*, -00 = 7r, the system state is initial- 
ized over a stable stationary point, and the correspond- 
ing orbit turns out to be trapped for P$ smaller than the 
maximum resonance width. 

These analytical predictions are in agreement with the 
FMA map shown in Fig. [6] Indeed, the bright green 
round curve in Fig. [6] corresponds to the separatrix of 
the resonance oj r = 2w0, which can also be visualized by 
plotting the average P$(t) for different anti-atom trajec- 
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field strength can be written as: 



1000.2 1000.4 1000.6 1000.8 
Time, s 



1001.2 1001.4 1001.6 



FIG. 8. Survival fraction f(t) = n(t) / n(t s hut) in the 
quadrupole trap under the assumption that T = 1 (solid) or 
J- — 200 (dashed). Anti-atoms are initialized with a \[E dE 
distribution and E < 550 mK (black) and with a flat dis- 
tribution and 350 mK < E < 550 mK (gray). The vertical 
line corresponds to a time £r = 1000.16 s at which the radial 
trapping potential drops down to (7 m m = 350 mK. The inset 
shows a zoomed in region for t < £r,. 



tories (Fig. [7| . If the anti-atom orbit is trapped within 
this resonance, (P^) = 0, while for anti-atoms outside 
of the resonance island, (P^) is finite. By crossing the 
separatrix, one would therefore expect to see a jump in 
(Pcf,). The analytical predictions for the resonance width 
and for shown in Fig. [7] for if>^ = and = tt, would 
seem to agree with the (P^) jumping near the actual sep- 
aratrix. 



3. Radial barrier shutdown 

In the previous section, based on the numerical simu- 
lation of dynamics governed by the frozen Hamiltonian, 
we inferred that, for our numerical example, a signifi- 
cant fraction of trapped anti-atoms with E > U m - ln have 
regular trajectories. This makes the assumption of tra- 
jectory ergodicity unjustified. However, it is still possible 
that the anti-atom gravitational mass can influence how 
anti-atoms escape as the radial potential well lowers with 
the decrease of the quadrupole coil current I q . Suppose 
that the shutdown of the quadrupole coil starts at t = to- 
If, for a fixed profile I q (t), the fraction of trapped anti- 
atoms f(t) = n(t)/n(t ) is different for different values of 
J 7 , one can use an experimental measurement of f(t) to 
infer bounds on the gravitational mass. In the following, 
we compare simulations of f(t) for T = 1 and J- = 200. 

We next identify the multipole field with the field 
B q (r,t) created by the quadrupole coils. Introducing 
a[t) = I q (t)/I q (t ) so that B q (r,t) = a(t) B q (r,t ), the 



\B\ 



Bl + Bl + B'l = \[B bz + B pz + aB qz {r,t )f 
[B pr + aB qr (r,t )} 2 + a 2 B 2 q<f> {r, t )| 



1/2 



B\ + 2aG + a 2 B2{r,t ), (21) 



where B\ = 



{Bbz + B pz ) 

{Bbz + B, 

radial and azimuthal components of vectors, respectively. 



- , pz)B qz (r,to), and the subscripts r 



pr , G — B pr B qr (r,t ) + 
denote the 



Equation (21) was implemented numerically by tabu- 
lating zeroth-order and second-order azimuthal harmon- 
ics of B\, G and B 2 (r,t ) independently. The simu- 
lated antihydrogen ensemble contained 64,000 anti-atoms 
with an energy distribution M{E) dE scaling like \[E dE 
[21 [211 H21 • All anti-atoms were initialized with energy 
below 550 mK, because any anti-atom with higher en- 
ergy leaves the device within 10 ms. We compared the 
loss of anti-atoms due to the quadrupole coil shutdown 
for T = 1 and T = 200. For the first t = 1000 sec- 
onds, the quadrupole coil is energized a(t < t Q ) — 1. 
Then, the quadrupole coil is turned off with a character- 
istic time scale on the order of one second. A choice of 
a(t) = exp[-2(i - t a ) 2 /(t -t + 0.8 s)] for t > t , similar 
to the reconstructed radial trapping potential shown in 
Fig. 3b of Ref. E 

The time-dependence of the fraction of anti-atoms re- 
maining trapped after the initiation of shutdown is shown 
in Fig. [8] According to Fig. [8j the dependencies f(t) 
calculated for T = 1 and T = 200 are virtually iden- 
tical. Introducing the moment of time « 1000.16 s 
when the radial potential barrier at z — drops down 
to t/mm, one observes that, while approximately 3% of 
anti-atoms escape the device prior to £r in a system with 
T = 1, about 2.5% of anti-atoms escape over the same 
time interval when T = 200. Note that if the anti-atom 
motion were ergodic, no anti-atom de-trapping would be 
observed until t = tn for the case where J- = 200. 

Of course, the fact that in ATRAP, about 10% of all 
annihilation events were detected before t = ip, [5] could 
be attributed to the fact that the actual anti-atom dis- 
tribution function might differ significantly from y/E dE. 
Additional simulations performed with a flat distribution 
function, containing only anti-atoms with energies in a 
range 350 mK < E < 550 mK, were shown to be in close 
agreement with results obtained for a M{E) dE oc y/E dE 
distribution and, in this case, the fraction of anti-atoms 
escaping before t = % reached 7%. A small deviation of 
0.5% between the graphs of f(t) for t < £r shown in Fig. [8] 
could, in principle, be detected in an experiment. Note, 
however, that we infer that only approximately 4 anti- 
hydrogen annihilations (with 5 expected cosmic events) 
were observed in total in Ref. [5] in the relevant time re- 
gion between t = to = 1000 s and i = t R = 1000.16 s. 
This count rate is at least two orders of magnitude lower 
than that necessary to resolve the differences between the 
curves shown in Fig. [8j 
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FIG. 9. Dependence of {z)/L on r for 100 mK antihydrogen 
atoms trapped in a vertical octupole trap with parameters 
similar to those of ALPHA (Sec. Ill B l with: (a) artificial 
separable potential U(r, z) = U (r, 070J+ (7(0, 0, z) - (7(0, 0, 0) 
(green solid) and (b) realistic trapping potential U (r, <f>, z) 
(blue solid). Two fitted exponential functions exp(— 2.6t — 
0.2) - 1 (green dotted) and 0.5(exp(-4.1r) - 1) (blue dotted) 
are shown for reference. 



We infer from these simulations that one cannot es- 
tablish a limit of J 7 < 200 using the technique described 
above with a vertical quadrupole trap. Note that the two 
distributions studied here are very different, but lead to 
the same conclusion. 

In the following section we turn our attention an im- 
proved technique. 



4- Axial barrier shutdown 

A natural alternate approach for measuring the anti- 
hydrogen gravitational mass involves lowering the axial 
trapping barrier in a vertical trap. Recall that if M > 0, 
the gravitational potential Mgz lowers the trapping po- 
tential at the bottom of the trap and raises it at the 
top, relative to the trap center; for M < 0, the trap- 
ping potential is lowered at the top and raised at the 
bottom. Assuming that currents in both coils are very 
nearly equal at each moment of time, and that the mag- 
netic field they produce is decreasing in magnitude suf- 
ficiently slowly, nearly all anti-atoms with M > will 
be expected to exit at the bottom of the trap, where the 
trapping potential is slightly lower (Fig. [IJ. For M neg- 
ative, antihydrogen would instead preferentially exit the 
trap at the top. Observing the vertical location of anti- 
hydrogen annihilations during slow shutdown of mirror 
coils may, therefore, be a useful experimental technique 
for quickly assessing the sign of M. Some preliminary 
estimates of the required shutdown time and a numerical 
simulation of such an experiment are discussed below. 



The characteristic adiabatic time-scale t*, on which 
the trapping potential should be lowered in order to de- 
termine the sign of M can, in principle, be estimated 
by analyzing the axial motion under the Hamiltonian 
Hn (jp z , z, t). Suppose that M > and consider an anti- 
atom which is about to cross the inner separatrix S\ pass- 
ing through the saddle point of the lower potential barrier 
(Fig. fll). Let T{Hn,t) be the period of the antihydrogen 
trajectory calculated for a frozen potential profile U(z), 
and let T be the smallest period of all orbits between 
two separatrices. If T is sufficiently large, the anti-atom 
may cross another separatrix S2 passing through the sad- 
dle point of the upper potential barrier, after the time 
At = (2MgL/U)r, where U is the trap depth and r is 
the actual field shutdown time. As a result, the proba- 
bility for such an anti-atom to leave the device at the top 
(p z > 0) will be approximately equal to the probability 
of leaving at the bottom (p z < 0) . On the other hand, if 
T -C At, nearly all antihydrogen crossing Si will leave 
the device at the bottom before reaching <S 2 - The field 
shutdown is then adiabatic if it occurs on a time-scale 
much larger than defined by: 



2MgL 
~1T~' 



T. 



(22) 



The value of T can be estimated by recalling that T{Hn ) 
goes to infinity (logarithmically) near both separatrices, 
and the minimum of T is therefore comparable to 



In 



4C7 
MgL 



2L 

77 



(23) 



where l z is the characteristic scale length of the axial 
confining potential. Assuming that M = m, this es- 
timate suggests that for 100 mK antihydrogen with a 
Gaussian distribution trapped in a device similar to AL- 
PHA, r* is expected to be of an order of the second. 
To verify this conjecture, numerical simulations of an- 
tihydrogen escaping from a trap with separable poten- 
tial U[r,z) = U r {r) + Uz(z), where U r (r) = J7(r,0,0), 
U z {z) = 17(0,0, s) - £7(0,0,0), and U{r,<p,z) = \iB 
is the confining potential of an ALPHA-like appara- 
tus described in Sec. |III B| were performed. Lowering 
the current in the mirror coils according to I m (t) — 



ImO e 



-(t-t )/T 



■ calculated the z coordinates of all sim- 



ulated annihilation events and plotted their average (z) 
as a function of r. As expected, this average annihilation 
position (z)(t), shown in Fig.|9j converges to the bottom 
of the trap — L as r goes to infinity. The characteristic 
time-scale of this dependence is on the order of a second, 
in agreement with the prediction for . 

If implemented, this experimental technique could po- 
tentially allow one to distinguish between T > a(r) and 
J- < —<x(t) with a(r) — > as r — > 00. Choosing a 
sufficiently large t, it might even allow one to distin- 
guish T > 1 from T < —1 for even 300 mK antihydro- 
gen atoms. Unfortunately, however, this proposed tech- 
nique would be very sensitive to possible deviations of 
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FIG. 10. Reverse cumulative averages (z)J of all 3.2- 10 5 an- 
nihilation events and corresponding confidence regions (blue 
for T = 300 mK, red for T = 100 mK and green for T = 
30 mK) obtained numerically for anti-atoms with T = 1 (bot- 
tom) and F = — 1 (top) in a vertical trap with r = 300 ms. 
The assumed particle distribution is \/i5exp(— E/kT). 
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FIG. 11. Reverse cumulative averages (z)^ of all 3.2- 10 5 an- 
nihilation events and corresponding confidence regions (blue 
for t — 50 ms, red for r = 100 ms and green for r = 300 ms) 
obtained numerically for 100 mK anti-atoms with F = 1 (bot- 
tom) and F = — 1 (top) in a vertical trap. The assumed 
particle distribution is \/Eexp(— E/kT). 



the actual trapping potential U(r, 4>, z) from its separa- 
ble approximation U(r,z). One consequence of the non- 
separability of U (r, <fr, z) is the emergence of stochastic 
layers near the separatrices. If the layers overlap, anti- 
atom dynamics within the region confined by S\ and S2 
will be stochastic, and the approximate expression de- 
rived for T will no longer be valid. On the other hand, 



a small non-separable field component SU = U — U may 
perturb low-energy antihydrogen trajectories as the mir- 
ror coils are being shut down. Indeed, if the oscillations 
of I z due to the perturbation SU are sufficiently strong 
and exceed the distance between the separatrices, anti- 
atoms will cross both of them numerous times. As a 
result, anti- atoms with M > will retain a finite prob- 
ability of leaving the device at the top, even if the field 
shutdown is infinitely slow. 

This effect can be observed by simulating 100 mK an- 
tihydrogen escape from a device with a realistic trap- 
ping potential U(r,<j),z). Now (z)(t) does not converge 
to — L as t — > 00, but instead becomes saturated at 
(z) rj —7.5 cm (Fig. [9]). Therefore, increasing the shut- 
down time r (beyond about 0.5 s in our case) does not 
necessarily lead to a substantial decrease of a(r) nor to 
improvement of the antihydrogen mass measurement. 

5. Accuracy of the gravitational mass measurement 

As discussed in the previous section, (z) of escaping 
particles measured in a vertical trap with de-energized 
mirror coils can be very sensitive to the gravitational 
mass M of the cold antihydrogen. This effect could, in 
principle, be used to determine the value of M, or simply 
check whether M is greater or smaller than zero. Here, 
we determine an accuracy of such a hypothetical test by 
comparing (z) for F — — 1 and F = 1 assuming that there 
are only 500 annihilation observations. To accomplish 
this, we follow the procedure discussed in Ref. @1 namely 
we calculate the reverse cumulative averages (z)* k (t) for 
640 sets Sfc, each containing 500 simulated annihilation 
events. This reverse cumulative average is defined as the 
average z of events occurring after the time t, i.e., 

(*)m = ( E E • ( 24 ) 

\iGs k ,ti<t / \i£s k ,ti<t ) 

Likely statistical fluctuations of {{z)* k } can then be visu- 
alized by plotting a confidence region [z\(t), 22 (i)], cho- 
sen in such a way that the intervals (—00, z%) and (z%, 00) 
each contain only 5% of the values of {{z)* k }. 

The simulations were performed for a vertical trap with 
parameters similar to those of the ALPHA apparatus. 
The confidence regions for J- = ±1, but different values 
of the antihydrogen temperature and shutdown times are 
shown in Figs. [T0| and [XT] In these simulations, the time 
profile of the current in the mirror coils was chosen to 
be I m (t) = I m oe~ t/T with t = 0.05 s, 0.1s, and 0.3s 
(we have chosen to = for simplicity). According to 
our results, a measurement of (z) for r = 50 ms and a 
\/Eexp(—E/kT) dE particle distribution with T up to 
at least 600 mK can be used to distinguish between F = 
— 1 and F = +1 hypotheses with a 95% confidence. In 
other words, a measurement (z) > (< 0) is inconsistent 
with F = +1 (J- = — 1) hypothesis since this average 
lies outside of the 95% (z) -confidence region simulated 
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for T = +1 (J 7 = —1). Note that calculating (z) of 
late-time events can further improve the accuracy of the 



method (Fig. 10). Simulations performed for a horizontal 
ALPHA trap suggest that a similar test on the sign of M 
can be accomplished only for cold plasmas and a large 
octupole coil shutdown time. Fixing the antihydrogen 
temperature at 30 mK, we see that two 95% confidence 
regions for T — ±1 intersect when r < 0.2 s. 



IV. CONCLUSIONS 

Measuring the ratio of the gravitational to inertial 
mass in neutral antihydrogen is possible in vertical and 
horizontal traps, but will require detailed simulations of 
the nonlinear dynamics of trapped anti-atoms, as these 
dynamics affect the nature of any signal of the gravita- 
tional interaction, and limit the accuracy with which it 
might be extracted. Our study of a vertical quadrupole 
trap based on the ATRAP experiment shows that the 
claimed experimental sensitivity is not realized with an 



experimental methodology inferred from Ref. [5J Surpris- 
ingly, insufficient stochasticity can limit schemes to mea- 
sure the gravitational mass of antimatter. In particular, 
because of a lack of ergodicity, radial shutdown in a ver- 
tical trap does not appear to offer much sensitivity to 
J- = M/m. The coupling of axial and transverse motions 
and the related notion of stochasticity of typical trajec- 
tories in phase space plays especially important roles in 
other measurement techniques as well. In some cases, 
Arnold diffusion [ITJ UHl [23] and other consequences of 
stochasticity can limit the precision with which gravi- 
tational interactions can be inferred. Systematic effects 
from small field errors and detector misalignments also 
need to be carefully understood. 
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